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Abstract. In high-energy nuclear collisions the degree of thermalization at the 
partonic level is a key issue. Due to their large mass, heavy-quarks and their 
participation in the collective flow of the QCD medium constitute a powerful 
tool to probe thermalization. We propose measuring azimuthal correlations of 
heavy-quark hadrons and products from their semi-leptonic decay. Modifications 
or even the complete absence of initially, e.g. in pp collisions, existing azimuthal 
correlations in Pb+Pb collisions might indicate thermalization at the partonic 
level. We present studies with PYTHIA for pp collisions at the top LHC energy 
using the two-particle transverse momentum correlator (Apt^i Apt^2) as a sensitive 
measure of azimuthal correlations. 



1. Introduction 



High-energy nuclear collisions offer the unique opportunity to probe highly excited 
nuclear matter in the laboratory. At sufficiently high temperature and/or energy 
density hadrons dissolve and quarks and gluons carrying color charge are the relevant 
degrees of freedom, commonly called a Quark Gluon Plasma (QGP). An essential 
difference between collisions of elementary particles on the one hand and heavy nuclei 
on the other hand is the development of collectivity in the latter. Collective flow 
of hadrons, especially the multi-strange hadrons (j) ^.nd fl, has been experimentally 
observed at RHIC [T] suggesting that collectivity dominantly develops in the early 
partonic stage, i.e. among quarks and gluons. Presently, the degree of thermalization 
among partons remains a crucial issue. 
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Figure 1. (Color online) 
Distribution of DD pairs as 
a function of relative az- 
imuth A(p from pp collisions 
at y/s = 14 TeV as calcu- 
lated by PYTHIA (v. 6.406) 
integrated over rapidity from 
flavor creation (red line) , 
flavor excitation (blue line), 
gluon splitting (green line), 
and their sum (black line). 
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Heavy-quark hadrons and their observables are of particular interest when 
addressing thermalization [5]. In a QGP, where chiral symmetry might be partially 
restored and consequently light-quarks (up, down, strange) obtain their small bare 
masses, heavy-quarks remain massive [3j. Heavy-quarks are mostly created in the 
early stage of the collision, since their mass is much larger than the maximum possible 
temperature of the QCD medium created while annihilation in the QCD medium is 
negligible |4] . Thus, heavy-quarks probe the entire history of a high-energy nuclear 
collision. Heavy-quarks participate in collective motion provided their interactions at 
the partonic level occur frequently. Thus, collective motion of heavy-quark hadrons 
is a powerful tool when addressing early thermalization of light quarks in high-energy 
nuclear collisions. 

1.1. DD angular correlations 

In strong interactions heavy-quarks are always created together with their anti-quark 
and are thus correlated. In collisions of elementary particles, these correlations survive 
the fragmentation process into hadrons to a large extent [3] and are thus observable 
e.g. in the angular distributions of pairs of D~ and D-mesons [5]. 

In high-energy coUisons of heavy nuclei, frequent interactions among partons 
(quarks and gluons) of the QCD medium and heavy-quarks would lead to a weakening 
and finally complete vanishing of these correlations. On the other hand, hadronic 
interactions at the late stage are insufficient to alter azimuthal correlation of Z? im- 
pairs [3]. It is frequent interactions distributing and randomizing the available (kinetic) 
energy and finally driving the system, i.e. light-quarks and gluons, to local thermal 
equilibrium. Thus, the decreasing magnitude when compared to pp collisions of 
heavy-quark correlations in high-energy coUisons of heavy nuclei would indicate early 
thermalization among partons. 

We employed the Monte Carlo event generator PYTHIA which reproduces 
experimentally observed correlations of D mesons at fixed target energies [5]. 
Parameters in PYTHIA were tuned to reproduce the NLO predictions O |H]. Our 
results on the yield of DD-pairs versus their relative azimuth for pp collisions at 
the top LHC energy of ^/s = 14 TeV are shown in Fig. [l] Our calculations at 
leading order (LO) contain flavor creation processes {qq — > QQ, gg — > QQ) and 
lead to an enhancement around w 180° , i.e. back-to-back . However, next-to- 
leading order (NLO) contributions such as flavor excitation {qQ qQ,gQ — > gQ) 
and gluon splitting {g QQ) become dominant at LHC energies and do not show 
pronounced correlations in the representation chosen in Fig. [T] leading to a rather weak 
dependence on relative azimuth. By choosing various momentum cuts, a somewhat 
stronger correlation was extracted [9]. 

In this paper, we introduce the two-particle transverse momentum correlator 
as a sensitive measure of heavy-quark correlations. This method has the following 
advantages: 

(i) the correlator is sensitive to non-statistical fluctuations, thus carving out any 
physical correlations and eliminating the need for various momentum-cuts in the 
analysis. This is essential when comparing results from different collision systems, 
e.g. pp and Pb+Pb. 

(ii) in case of physically uncorrelated candidate-pairs (e.g. background), the extracted 
value for the correlator is null, thus providing a reliable baseline. 
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2. Employing the two-particle transverse momentum correlator 

The occurence of non-statistical fluctuations of the event-by-event mean transverse 
momentum Mpt goes along with correlations among the transverse momenta of particle 
pairs. Such correlations were successfully extracted from experimental data employing 
the two-particle transverse momentum correlator [101 [H]. For D- and Z?-mcsons, this 
leads to 

(ApM,Ap,.)(--)=^5J^ (1) 
where is the pt convariance: 

j=i 

where pu and ptj are the transverse momentum pt for i*'* and j*'' particle of an event 
of D- and D-mesons respectively, pt is the inclusive mean transverse momentum 
averaged over all particles of all events of D and D, X]fc=i A^l"^"^ the total number of 
DD pairs and n^y the total number of pp collisions. 




Figure 2. (Color online) Distribution of the momentum correlator (Apt i, Apj 2) 
of DD pairs as a function of relative azimuth A0 at mid-rapidity (circles), 
integrated over rapidity (squares) and for background using the mixed event 
method (triangles) from pp collisions at y/s = 14 TeV as calculated by PYTHIA 
(v. 6.406). The lines are drawn to guide the eye. 

We studied the scale dependence of pt correlations in azimuthal space by 
calculating the correlator in bins of the relative azimuthal separation Acj) of particle 
pairs. For the case of independent particle emission from a single parent distribution, 
the correlator {Apt^i, Apt^2) vanishes. 

The DD momentum correlator {Apt.i, Apt^2) as a function of relative azimuth 
Acj) is shown in Fig. [2]for pp collisions at y/s = 14 TeV exhibiting a rich structure. We 
observe an enhancement at small azimuth from gluon splitting processes, while flavor 
creation of cc-quark pairs leads to a pronounced distribution at backward angels. 
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Figure 3. (Color online) Relative contribution to the production of D- and D- 
mesons (left) and B— and B— mesons (right) as a function of transverse momentum 
for the processes of flavor creation (red line) , flavor excitation (blue line) and 
gluon splitting (green line) from p — p collisions at y's = 14 TeV as calculated by 
PYTHIA (v. 6.406). 



Flavor excitation processes involve a larger number of gluons leading to a rather flat 
distribution At mid-rapidity, the correlations are even stronger, reflecting the harder 
Pt spectrum of DD meson-pairs when compared to the full rapidity range. Integrating 
the correlator over all azimuth, we get {Apts,^Pt,2) = 0.199 ± 0.006 GeV^/c^ 
which corresponds to the normalized dynamical fluctuation Spt [12 of ~ 30% in 
p^. This large value implies a strong correlation when compared to ~ 1% observed for 
unidentified charged particles in central collisions at SPS and RHIC [121 [El [14] . To 
mimic combinatorial background which is always present in the experiment, we applied 
the correlator to D- and D- mesons from different pp collisions, which are physically 
uncorrelated. This results in a value of (Apf.i, Apt^2) consistent with zero. Therefore 
the correlator allows for a clear distinction between the case were correlations are 
present (different from zero) or absent (equal to zero) in contrast to the method 
described in [3 . The different QCD processes contributing to heavy-quark productions 
have strong energy dependences. Figure |3] shows the relative contribution to the 
production of D~ and D-mesons (left) and B- and S-mesons (right) as a function of 
transverse momentum. The production of Z?-mesons is dominated by gluon splitting 
and flavor excitation processes while the contribution from flavor creation is about 10% 
at low momentum and increases up to 20% at larger momentum. At RHIC energies, 
the contribution from gluon splitting was estimated by the STAR experiment to be 
less than 10% (TS]. On the other hand, the production of _B-mesons is dominated 
by flavor creation and flavor excitation with a small contribution from gluon splitting 
below 10% and an overall weak dependence on transverse momentum. 

As shown above, the initial correlations of cc-quark pairs survive the 
fragmentation process to a large extent. However experimentally, full kinematic 
reconstruction of D-mesons from topological decays suffer from small efficiencies 
resulting in low statistics, especially when pairs of D-mesons are considered where the 
reconstruction efficiency enters quadratically. To circumvent this fact, we considered 
electrons (positrons) from semi-leptonic decays of D- and B-mesons with a branching 
ratio to electrons of 10% and 11%, respectively. 

Figure|4]shows the momentum correlator {Apt^i, Apt^2) of e"'"e~-pairs from decays 
of DD (squares) and BB-mesons pairs (triangles) as a function of their relative 
azimuth integrated over the full rapidity range. We chose a momentum cut 
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Figure 4. (Color online) Distribution of the momentum correlator 
(Apt J, Apt, 2) of e+e~ pairs stemming from semi-leptonic decays of pairs of SB- 
mesons (squares) and DD-mesons (triangles) as a function of relative azimuth 
Alp at mid-rapidity in pp collisions at y/s = 14 TeV as calculated by PYTHIA (v. 
6.406). The lines are drawn to guide the eye. 

of pt > IGeV/c allowing for clean electron identification at reasonable background 
and still sufficient statistics in the ALICE experiment using the Transition Radiation 
Detector. We observe a correlation at forward angles due to gluon splitting and a more 
pronounced correlation at away-side angles from flavor creation processes indicating 
that the initial correlations among heavy-quarks and their anti-quark even survive 
the semi-leptonic decay into electrons (positrons) making an experimental observation 
at the LHC possible. We note that the correlation is stronger for electrons from semi- 
leptonic decays of _B-mesons compared to D-mesons. This is in agreement with results 
from the STAR experiment at RHIC [16]. 

When integrating over full azimuth for D-mesons, we get a value for the correlator 
of {Apt^i, Apt^2) — 0.007 ± 0.001 GeV^/c^, which corresponds to the normalized 
dynamical fluctuation Spt of ^ 12%. Hence, more than half of the correlation strength 
observed for fully reconstructed Z?-meson pairs is still present in pairs of electrons and 
positrons from semi-leptonic decays. 

3. Conclusions and outlook 

In summary, we propose measuring correlations of heavy-quarks and their 
modifications in Pb-|-Pb collisions compared to pp collisions as a sensitive probe of 
thermalization at the early partonic stage. We applied the momentum correlator as 
a precise and normalized measure. Different QCD processes contributing to heavy- 
quark production are identified in the azimuthal distribution of the correlator. Our 
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results indicate that these correlations born at the creation of pairs of heavy-quarks 
and their anti-quarks survive the fragmentation in pp collisions into hadrons and even 
semi-leptonic decays into electrons to a large extent. Thus, experimental observation 
of these correlations in pp and Pb-|-Pb collisions at the LHC seems feasible. 
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